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ABSTRACT 

We study the width of the semiclassical 
phase distribution of a quantum state in its 
dependence on the average number of photons 
<m> in this state. As a measure of phase 
noise, we choose the width Atp of the best 
Gaussian approximation to the dominant 
peak of this probability curve. For a coherent 
state this width decreases with the square root 
of (m) , whereas for a truncated phase state it 
decreases linearly with increasing (m). For 
an optimal phase state, Ay decreases 
exponentially — but so does the area “caught” 
underneath the peak: All the probability is 

stored in the broad wings of the distribution. 

I. INTRODUCTION 

The ultimate quantum limit in the goal of 
optically detecting gravitational waves is to 
operate a Michelson interferometer with light 
in a quantum state that minimizes the phase 
noise at a given mean number of photons 
(Ref. 1). But what is a good measure for phase 
noise? Should we consider the inverse peak 
height of the probability distribution (Ref. 2) 
— the so-called reciprocal likelihood — or 
perhaps the second moment of the phase 
distribution (Ref. 3) or even the periodic 
measure advocated in Refs. 4 and 5? These 
are all based on the idea of a phase 
distribution — but we recall that this in itself 
is not a trivial construction since the concept 
of a Hermitian phase operator is not without 
complications (Ref. 6). 


In the present paper we therefore start 
from the semiclassical phase distribution 
W(p[ I y)] of a quantum state I y) (Refs. 7 and 
8). We consider states that have a single 
pronounced maximum of W(p at a phase value 
that, without loss of generality, we take at (p = 
0. We approximate (Ref. 9) this peak by a 
Gaussian distribution with an identical 
height of W<p=o. The distribution's width Ay 
is determined by the curvature of W<p at <p=0, 
together with W<p=o. The examples of a 
coherent state I y C oh) and a truncated phase 
state ly p ) will illustrate this scheme: Their 
widths Ay [ I y C oh>] and A(p [ I y p >] decrease as 
the square root and linearly with increasing 
average photon number (m), respectively. In 
the case of the optimal phase state — the state 
(Ref. 2) that minimizes the reciprocal 
likelihood — we find (Ref. 9) even an 
exponential decay for Acp. However, in 
contrast to the coherent or the truncated phase 
state, the area underneath this maximum is 
not normalized but decays as well (Ref. 10). 
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II. FROM PHASE FUNCTIONAL TO 
GAUSSIAN-APPROXIMATED 
PHASE DISTRIBUTION 
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In this section we start from the 
semiclassical phase amplitude functional 
w<p [l\|/)l of a quantum state I y) and derive a 
Gaussian approximation to the dominant 
maximum of this probability curve. 

In the semiclassical limit the phase 
distribution [ I \jr)] of a quantum state 
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Here primes denote differentiation with 
respect to <p and we have used the property 
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represented as a superposition of photon 
number states Im), with expansion 
coefficients \\f = (m I vj/) and the 

normalization ^ such that (\\f\ y) = 1, follows 
from the phase functional 
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For the sake of simplicity, we consider in this 
article only quantum states such that = 
* 

> 0. Hence, for the phase value <p = 0 the 

terms in the sum of Eq. (2.2b) add 
constructively, whereas for cp * 0 
cancellations occur. This results in a 
maximum at <p = 0. 

An approximate analytical expression for 
W(p [l\j/>] in the neighborhood of the origin 
follows from an expansion of W<p into a 
Taylor series around <p = 0, that is 


following from Eq. (2.2b). With the help of 
Eq. (2.3), we arrive at 
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where the width A<p of this Gaussian is given 
by 
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We emphasize that this procedure is valid for 
any state whose phase probability Wq> [I y)] 
enjoys a maximum at cp = 0. 

The area underneath this Gaussian- 
approximated peak reads as 
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cp=0 

This clearly demonstrates that the Gaussian 
fit of the peak, Eq. (2.4), is different from the 
properly normalized Gaussian G((p) = ti -1 ' 2 

o (peak) 

(Acp) -1 exp [-(cp/Acp)^]. Whereas is 

tailored to have a height identical to 
W(p = o [ I V)1 , the height of G, that is, rc -1 ' 2 
(Acp) -1 , adjusts itself to the width of the 
Gaussian, so as to keep the area normalized. 


recently proposed in Ref. 2. 

These states are normalized to unity. 
For a coherent state, Eq. (3.1), we arrive at 



III. PHASE NOISE AND AVERAGE 
NUMBER OF PHOTONS 

In this section we apply the Gaussian 
approximation, Eq. (2.4a), in order to discuss 
the width Acp , Eq. (2.4b) of the phase 
distribution of: 

(i) a coherent state of large displacement 
a » 1 , 



(3.1) 


(ii) a truncated phase state 


m 

/ \-l/2 «-? 

| ¥ )=(m 0+ lj (3 ‘ 2) 
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and 


(iii) the optimal phase state 


= 1 (3.4) 


where we have replaced the summation by an 
integration. For the truncated phase state 
I \|/p) we find directly 

-1 m o 

<yy = ( m 0 +1 ) X 1 - 

m=0 

whereas for the optimal state Eq. (3.3), 

the normalization constant is given 

implicitly by 

m o 

(Vg'Vg) = ^ X (1+m) 



that is, (Ref. 11) 
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and the width Acp , given by Eq. (2.4b) reduces 
to 


= $ 
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Here, 


C(S,D) = £(D + k) S 
k=0 


denotes the generalized Riemann zeta 
function (Ref. 12). 


A* Coherent State 

For the coherent state, Eq. (3.1), the 
expansion coefficients read 
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Here, we have once again replaced 
summations over m by integrations. We 
evaluate the Gaussian integrals most 
economically by applying the symmetry of 
\| / with respect to or before performing the 

integrals. This yields 
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This result reduces Eq. (3.6) to 
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We now evaluate the area underneath 
the Gaussian approximation, Eq. (2.5), for 
the coherent state, Eq. (3.1). The maximum 
value of the phase distribution reads 




(3.7a) 


The average number of photons (m) 
follows from the normalization condition, 
Eq. (3.4), as 


(m> 



= Jdm(m-a 2 )<n 
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and hence 


and hence the area of the Gaussian is 
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Thus the Gaussian approximation, Eq. (2.4) 
for a coherent state is properly normalized. 
Its width Acp [ I Vcoh)] decreases linearly with 
a but its height W<p=o [ I V co h>^ increases 

linearly with a ~ V (m) to leave the area 
Aoauss constant. 


B. Truncated Phase State 

We now turn to the discussion of the 

truncated phase state I \j/ p ) , Eq. (3.2). 

According to this equation, the expansion 

coefficients vp read 
T m 


Acp 




(3.7b) 


{ 1 , for 0 < m < 
0 , for m > m 

o 


for large <m) . Asymptotically, the width 
Acp [ I Y C oh^ P^ ase distribution of a 

coherent state decreases inversely as the 
square root of the average photon number. 


which reduces Acp [ I \p p )] to 
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we find 
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again for large (m) . The width Acp [ I y p >] of 

the phase distribution of the truncated phase 
state decreases linearly with the average 
photon number (m) . 

The value of Wq> at (p=0 reads 
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Here, we have made use of the summation “ + l 

formulae (Ref. 11 ) 
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When we express m 0 in terms of the 
average number of photons 


= - (m> , (3.13) 

71 

and yields for the area underneath the 
Gaussian approximation, Eq. ( 2 . 4 ), of the 
truncated phase state, Eq. (3.2), 


AGausst ' V]=^^<p[ ' V] * \ =Q [ 1 V] 



i.e., almost perfect normalization. Again, we 
note from Eqs. (3.10) and (3.13) that the width 
A 9 decreases with m 0 in the same manner 
as the maximum height W<p = o increases — 
keeping the area Acauss normalized and, 
more importantly, keeping it independent of 
m 0 . 
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C. Optimal Phase State 


When we use the summation formulae, (Ref. 

ID 


We conclude our illustration of the 
Gaussian-approximated phase distribution by 
applying it to the optimal phase state I v 3 > , 

Eq. (3.3), which enjoys the expansion 
coefficients 
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In the large m Q limit we arrive at 


(3.15) 
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and with Eq. (3.15) we then have 


where we have defined y = e c in a standard 
notation. In determining the remainder, we 
have applied in the last step of Eq. (3.17a) the 
asymptotic expression for ?£ , Eq. (3.5). With 
the help of Eqs. (3.17a) and (3.17b), Eq. (3.16) 
reads 
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The average number of photons in this state 
reads 
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which shows that the Gaussian-approximated 
width A<p of the phase state I \\f 3 ) decreases 
exponentially with the average number of 
photons. The height 
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of the distribution increases only 
quadratically . Here, we have also made use 
of Eq. (3.17a). Hence, the phase probability or 
area caught underneath the peak, 
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decreases rapidly to zero with increasing 
average photon number. All the probability is 
stored in the broad wings of the distribution. 


IV, SUMMARY AND OUTLOOK: A NEW 
VARIATIONAL PROBLEM 


In this article we have presented a 
Gaussian approximation to the maximum of 
the semiclassical phase distribution of an 
arbitrary quantum state. We have illustrated 
this scheme using the example of a coherent 
state, a truncated phase state, and the 
intriguing optimal phase state. Our main 
results are summarized in Table 1. The 
Gaussian-approximated width Acp [ I Vcoh)L Eq. 
(3.7), of the coherent state I Vcoh) decreases 
as the square root of the average number of 
photons, whereas for a truncated phase state 
Vp > the width Acp [|y p >] , Eq. (3.12), decreases 
linearly with <m>. In both cases the Gaussian 
is properly normalized, that is, the probability 
caught underneath the peak is almost 


identical to unity, Eqs. (3.9) and (3.14), and 
independent of (m). 

The situation is quite different for the 
optimal phase state I Vs)- Here, Acp [ I Vs)L Eq. 
(3.18), decreases exponentially with (m), but 
the maximum W(p = 0 NVs)L Eq. (3.19), 
increases only quadratically with (m), 
leading to vanishing probability in the peak. 
All probability in this case is in the broad 
wings of the distribution, as is discussed in 
detail in Ref. 9. 

We conclude by noting that the 
Gaussian approximation might lead to 
insight into questions such as the existence of 
a lower bound of Acp for a given fixed 
number of photons <m>. For that purpose we 
would like all the probability to reside in the 
peak, that is, Aoauss = From Eq. (2.5), we 

find 


which, when substituted into Eq. (2.4b), yields 
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Two strategies offer themselves: (i) Use 

appropriate inequalities to rewrite the 
expression in square brackets in Eq. (4.1) in 
terms of the average number of photons and 
its variance. This might lead to a lower 
bound of Acp . (ii) Minimize Acp 3 , that is, 
maximize the expression in square brackets, 
Eq. (4.1), under the constraint of constant (m) 
and phase state normalized to unity. 
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Table 1. Gaussian approximation W<p = W<p = o exp [ — (cp/A(p)^] for dominant maximum of the phase 
distribution W 9 = ( 27 c)- 1 -H 2 | =0 exp (im<p) | 2 for a coherent state I y coh > , a 

truncated phase state I y p > and an optimal phase state I y s > • 
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